Structure of the superconducting state in a fully frustrated wire network 

with dice lattice geometry 



O 

o 

(N 

-t— > 
o 

o 

(N 
(N 

o 
o 

u 

Oh' 



O 

o 



> 
oo 
oo 
m 

o 

o 

-I— > 



I 

o 
o 



X 



S. E. Korshunov^'-^ and B. Dougot 

^L. D. Landau Institute for Theoretical Physics, Kosygina 2, Moscow 119334, Russia 
^Laboratoire de Physique Theorique et Hautes Energies, CNRS UMR 7589, 
Universite Paris VI and VII, 4 place Jussieu, 75252 Paris Cedex 05, France 
(Dated: October 22, 2004) 

The superconducting state in a fully frustrated wire network with the dice lattice geometry is 
investigated in the vicinity of the transition temperature. We express the projection of the Ginzburg- 
Landau free energy functional on its unstable subspace in terms of variables defined on the triangular 
sublattice of sixfold coordinated sites. For the resulting effective model, we construct a large class of 
degenerate equilibrium configurations, which are in one to one correspondence with ground states 
of the fully frustrated XY model with a dice lattice. The entropy of this set of states is proportional 
to the linear size of the system. Finally, we show that magnetic interactions between currents pro- 
vide a degeneracy lifting mechanism and find the structure of the periodic state selected by these 
interactions. 

PACS numbers: 74.81. Fa, 64.60.Cn, 05.20.-y 



I. INTRODUCTION 

The concept of frustration has been a common link 
among various problems in statistical mechanics for the 
past two decades at least. Even in the absence of dis- 
order, it often results in a phenomenon of competition 
between several degenerate ground-states. Superconduct- 
ing wire networks provide a very appealing class of sys- 
tems where many subtle effects induced by frustration 
can be observed experimentally and analyzed theoreti- 
cally [aaaii|6,i 8, 9]. 

For simple regular networks, a natural parameter char- 
acterizing the strength of the frustration is the ratio 
/ = <I>/$o where $ is the external magnetic flux through 
an elementary plaquette of the lattice and $o = /ic/2e is 
the superconducting flux quantum. For an ideal network 
of very thin wires, all physical properties are expected to 
be periodic functions of /, all integer values being equiv- 
alent. In this case, the maximal frustration is obtained 
when / reaches half-integer values. Such fluxes are in- 
teresting because already for a single loop, they provide 
two equivalent ground states, distinguished by the orien- 
tation of the supercurrent flowing around the loop. For 
more complex geometries, two adjacent loops (sharing a 
common link) have a lower free energy when the cur- 
rents in them flow in opposite directions. The possibility 
to fulfill this requirement for any pair of adjacent loops 
is a geometrical property of a given lattice, which allows 
one to be sure about the structure of the superconduct- 
ing state without any additional analysis. This clearly 
holds for a square lattice ^3 (where vortices of alternat- 
ing signs form a regular checkerboard pattern) or for a 
triangular lattice. 

In recent years, network geometries which do not sat- 
isfy this criterion have received a lot of attention. The 
most studied examples are the h oney con ib fllin^lTsLll^ 
m , the kamme M MM M IS 111 EaUjTand the 
dice miHilllallllla lattices. On the honeycomb 



lattice, the discrete degeneracy of the classical ground 
states in fully frustrated superconducting wire networks 
or Josephson junction arrays can be described in terms of 
formation of zero-energy domain walls in parallel to each 
other ^13j , the residual entropy of such system being pro- 
portional to its linear size. Experimentally, a cusp-like 
local maximum in the superconducting transition tem- 
perature Tc is observed as the external magnetic field is 
varied around the value corresponding to / = 1/2 jl^ . 
This behavior has been interpreted as an evidence for a 
degeneracy lifting mechanism which selects a commensu- 
rate ordered pattern of vortices |Q| . 

For the kagome lattice, the residual entropy of classical 
ground states is much larger, since it is proportional to 
the whole network area [Sg . The experimental situ- 
ation at /= 1/2 is not as clearcut as for the honeycomb 
lattice, since the shape of Tc versus magnetic field curves 
near / = 1/2 depends on the resistive criterion chosen to 
determine Tc [I3 , or on the superconducting metal (alu- 
minium versus niobium for instance 13). Theoretically, 
various degeneracy lifting mechanisms have been studied 
in detail by Park and Huse [20| . 

On the dice lattice (see Fig. P), the residual entropy 
is proportional to the system linear size |2^, as for the 
honeycomb lattice. Experimentally, magnetic decoration 
experiments 12^ have found a highly disordered vor- 
tex pattern, with a vortex correlation length comparable 
to the lattice spacing. Numerical simulations 27] of the 
corresponding XY model support the picture proposed 
in Ref. 26 for the ground states, but also demonstrate, at 
low temperatures, an unusually slow relaxation of energy, 
as well as aging of phase correlation functions. 

In this article we consider maximally frustrated super- 
conducting networks on a dice lattice, in the immediate 
vicinity of the superconducting transition temperature. 
In this limit, the amplitude of the superconducting order- 
parameter is not necessarily uniform, and it is appropri- 
ate to use a generalization of the approach introduced by 
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Abrikosov (23 in his first prediction of vortex lattices in 
type II superconductors. The main idea underlying this 
approximation is that at Tc{H), only a small fraction of 
the eigenmodes of the linearized Ginzburg-Landau equa- 
tions become unstable. Abrikosov developed a variational 
procedure where the superconducting order-parameter 
is constrained to remain in this unstable subspace (for 
T Tc{H) this procedure is asymptotically exact). Min- 
imizing the quartic term in the Ginzburg-Landau free 
energy functional yields then periodic vortex lattice so- 
lutions. 

Adopting this approach to a dice lattice network is 
quite interesting since the corresponding eigenmodes (for 
/ = 1/2) have unusually high degeneracy and exhibit the 
unexpected property of an extreme form of spacial local- 
ization. It is indeed possible to construct an eigenfunction 
basis for which each member is non-vanishing only on a 
finite cluster '55, 'SC/I . This phenomenon arises from the 
Aharonov-Bohm interference effect which is magnified in 
the geometry of the dice lattice, and these Aharonov- 
Bohm cages have been evidenced experimentally by the 
observation of magnetoresistance oscillations in ballistic 
semiconductor networks [3l| with the flux period hc/e 
per elementary loop. 

The main result of the present study is that in max- 
imally frustrated superconducting wire network non- 
linear effects select a class of order-parameter configu- 
rations in one to one correspondence with the ground 
states of the fully frustrated XY model with the same 
geometry j2£j| , which may be viewed as a low temperature 
approximation for the Ginzburg-Landau model ignoring 
amplitude variations (London limit). However, the inclu- 
sion into analysis of the magnetic energy leads to the re- 
moval of the accidental degeneracy and selection of one of 
the periodic states minimizing the Ginzburg-Landau free 
energy. The same state has the lowest free energy also at 
lower temperatures (down to London limit), as well as in 
Josephson junction arrays with the same geometry. 

In sections nil and HTll we express the Ginzburg-Landau 
functional for a fully frustrated dice lattice wire network 




FIG. 1: Finite cluster with the dice lattice geometry. 



after projection on the subspace of unstable modes in 
terms of complex variables defined on the triangular sub- 
lattice of sixfold coordinated sites. Section Hvl describes 
the construction of periodic equilibrium states for this 
effective problem, and their extension to a larger class 
of degenerate states whose precise connection with those 
proposed for the corresponding XY model is established. 
Finally, sectionlvl investigates a degeneracy lifting due to 
magnetic interaction between currents. 



II. HARMONIC CONTRIBUTION TO FREE 
ENERGY 

A. A single wire 

In the framework of the Ginzburg-Landau approxima- 
tion the free energy of a thin superconducting wire, F^^^, 
can be written as the sum of the two terms, 

dx {--|A(a;)p+ (1) 

2 



F, 



(2) 



and 



F, 



(4) 



. d 2ti 
dx <i>n 



^11 (x) 



A(x) 



dx\/^{x)\'^ , 



(2) 



describing, respectively, the harmonic and the fourth- 
order contributions to F^^^. Here a oc Tco — T, (3 and 
7 are the coefficients of the Ginzburg-Landau expansion, 
L is the length and Tcq the mean field transition tem- 
perature of the wire, A{x) is the superconducting order 
parameter as a function of the coordinate x along the 
wire, A||(a;) is the projection of the vector potential on 
the wire and $o — hc/2e is the superconducting flux 
quantum. 

At the point of phase transition |A(x)| 

(2) 

can be neglected in comparison with F^-^^.^. For the given 



and Fj-jl 



values of A(a;) at the ends of the wire, 

A(0) = Ao , A(i) = Ai , 
the minimum of F^^l^ is achieved when 



Aix) 



Ao sin ^ + Ai sin j exp(-iAoi; 



(3) 



exp ia{x) 



smrj 



where rj = L/£^, 



(4) 



(5) 



is 



lo the temperature dependent correlation length [here 
^ ~ ^(r = 0)], the function a{x) is defined by the inte- 
gral 



2tt r 

a(x) — -— / dx' A\\{x') , 



(6) 
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whereas ^oi is the value of this integral for the whole 
wire, = a{L). 

Substitution of Eq. Q into the expression for the su- 
perconducting current in the wire, 



lix) 



2e 



7 Re 



dx $0 



(7) 



shows that the value of the current is constant along the 
wire and is given by 



h sm rj '■ 



(8) 



On the other hand substitution of Eq. into Eq. 
gives a simple quadratic form of Ao and Ai [f^: 

FlZi^o,^i,Ao,) - F2.[cos77(|Aop + jAip) - (Ao|Ai) 



where F2 — 7/(2^ sin 77) and 



(9) 



(A,|Afe) = A,A 



+ A*Afce-*^^'= . (10) 



An arbitrary network 



The function F^j^ defined by Eq. @ can be then used 
to express the harmonic part of a free energy of a super- 

(2) 

conducting wire network Fnw in terms of the values of 
the superconducting order parameter Aj in its nodes j. 



i?(2) 



nw - X! -^whl (^i , Afc , Ajfe ) 



(11) 



Here the summation is performed over all links (jk) of a 
network. In the following we assume that all the links are 
identical and, therefore, the function F^^^.^{Aj , Ak, Ajk) 
is the same for all the links. 

In the case of a network formed by identical plaque- 
ttes it is convenient to express the value of perpendicular 
external magnetic field H in terms of the number of flux 
quanta per single plaquette: / = HS/^o (here S is the 
area of a plaquette). Then the directed summation of the 
variables Ajk = —Akj along the perimeter of each plaque- 
tte in positive direction (denoted below as ^q) should 
give 



jk 



2nf 



(12) 



From the form of Eq. ^ it is evident that the shift of 
/ by an integer or its reflection with respect to / = 1/2 
(/ — > 1 — /) do not change the form of the expression for 
free energy (or can be taken care of by a redefinition of 
variables), and, therefore, it is sufficient to analyze the 
interval < < 1/2. By the analogy with frustrated 
XY models [l^ a network with the maximal irreducible 
value of /, that is with / = 1/2, can be called a fully 
frustrated network. 



When fluctuations are completely neglected, the mag- 
netic field dependence of the superconducting transition 
temperature in a network Tc{f) can be found by looking 
when (with the decrease of temperature) the quadratic 
form defined by Eqs. lO-dJ looses its positiveness. To 

this end one has to analyze the system of equations ob- 

(2) 

tained by the variation of Fnw with respect to A 



^ [Afc cos 77 - 



A,e 



iAi 



(13) 



where j{k) are the nodes connected with k by the links 
of a network (in the following we call them the nearest 
neighbors of k) . The same equations can be derived P, Q 
directly in the framework of the continuous description 
without explicit calculation of -^^[^(Aj, A/j, Ajk). Multi- 
plication of Eq. (|13|l by A^, with subsequent extraction 
of the imaginary part allow to obtain the current conser- 
vation equation, 



3=j(k) 







(14) 



The form of Eq. (|13|l coincides with that of the 
Schrodinger equation for a single electron hopping be- 
tween the sites of the lattice with the same geometry in 
the presence of external magnetic field. As a consequence, 
Tc{f) can be related with eo{f), the largest eigenvalue of 
the Schrodinger equation in the same field. For a network 
whose nodes are all characterized by the same coordina- 
tion number z this relation can be written as 



Tco 



C eo(/) 

— arccos 

L z 



(15) 



Starting from the work of Hofstadter (who consid- 
ered the case of a square lattice), the spectrum of the 
Schrodinger equation for a single electron hopping prob- 
lem in the presence of external magnetic field has been 
extensively studied for various types of two-dimensional 
lattices includin g tr i ang ular |33| . honeycomb 34], dice 
[2^ and fcaoome l35l| lattices. 

The structure of the superconducting state in the net- 
work just below Tc{f) is determined by the structure of 
the eigenfunction corresponding to eo(/) @- The condi- 
tions for the applicability of the mean field approach are 
discussed in Appendix A. 



C. A network with a dice lattice geometry 

Dice lattice |3^ |33 is formed by the sites with the 
coordination numbers 3 and 6 in such a way that each 
bond connects two sites with different coordination num- 
bers (see Fig. 1). Below, when discussing a dice lattice, 
we denote the three-fold coordinated sites k and the six- 
fold coordinated sites j. Thus, the bond (jk) of a dice 
lattice connects the six-fold coordinated site j with the 
three- fold coordinated site k. 
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When considering the problem on a dice lattice it is 
convenient to simplify the quadratic form by mini- 
mizing it with respect to all variables defined on the 
three- fold coordinated sites. Substitution [from Eq. 
of 



. 1 T t 

3 cos r] ^-^ 



(16) 



into Eqs. (P)- ((TT|l then gives: 



f£) = ^V[(3cos^-l)(|A,f + |A,,| 



3 cos 77 

ui') 



2cos(7r/) (e'^^^' AjA*, -I- c.c 



•)] 



(17) 



where the summation is performed over the pairs {jj') of 
nearest neighbors on the triangular lattice formed by the 
six-fold coordinated sites, whereas variables 

= [{A,k' + Ak>f ) + {A,k" + A^ny )] /2 (18) 

(where fc' and k" are the two three-fold coordinated sites 
belonging to the same rhombus as j and j') are the av- 
erages of Ajji on the two shortest paths on a network 
connecting the nodes i and j'. It follows from Eq. H12|) 
that the variables A^y 
straint 



-Aj/j have to satisfy the con- 



(19) 



where N is the number of the six-fold coordinated sites in 
the network with appropriately chosen periodic boundary 
conditions (the total number of sites being SA*") . Accord- 
ingly, to find the structure of the superconducting state 
in a fully frustrated wire network with the dice lattice ge- 
ometry (which is the main subject of this work) one has 
to minimize the fourth-order contribution to free energy, 

^it^=E^wt(^-'^''-'^.-fc)' (23) 

[where A^ is given by Eq. H16() ] with respect to the whole 
set of the variables A^ satisfying the constraint H22|l . 

which fixes also the value of Fnw . For < 77 — 77c ^ 1 

« _127VF2A2(sin77,)(77 - r;,) , (24) 

where we have kept only the lowest order term of the 
expansion with respect to rj — rjc- 

At the conceptual level, this task can be considered as 
analogous to finding the structure of the vortex lattice 
which minimizes the forth-order contribution to the free 
energy of a bulk superconductor just below Hc2- In this 
problem (first analyzed by Abrikosov 29]), the harmonic 
contribution to free energy is degenerate with respect to 
a huge number of continuous variables, the positions of 
the order parameter singularities, whereas in the present 
problem a huge continuous degeneracy of the harmonic 
problem is related with variables A-,-. 



on all plaquettes of the triangular lattice. The form of 
Eqs. (jni) and ^ suggests that for < / < 1/2 the 
problem of finding Tc{f) on a dice lattice is reduced to 
analogous problem on a triangular lattice with / multi- 
plied by 3/2 and a different value of rj. Accordingly, the 
relation between the critical temperatures (expressed in 
terms of rj) in the two cases is given by 



3cos^77,(/) - 1 = 2cos(7r/)cos?7f (3//2) 



(20) 



Analogous relation between the single electron spectra 
on dice and triangular lattices has been derived in Ref. 

m 

Quite remarkably, for f = 1/2 the non-diagonal cou- 
pling in Eq. H17|l completely disappears, which allows im- 
mediately to conclude that 



ric (1/2) = arccos(l/\/3) w 0.9553 



(21) 



This absence of coupling between different variables A^ 
can be understood as a manifestation of the extremely 
localized nature of the highly degenerate eigenfunctions 
[2^ corresponding to the largest eigenvalue of the single 
electron Hamiltonian. 

(2) 

As a consequence, for f — 1/2 the value of FnJ turns 
out to be exactly the same for any set of the variables 
Aj satisfying the normalization condition 



(22) 



III. FOURTH-ORDER CONTRIBUTION TO 
FREE ENERGY 

A. A single wire 

Substitution of Eq. Q into Eq. Q describing the 
fourth-order contribution to the free energy of a super- 
conducting wire gives 



FlZi^„Ak,A,,) = F4[/4-(|A,|4-K|Afe|4) 

+ 2/3-(|A,|2 + |Afc|2)(A,|A 
+ /2- (2|A,|2|A,f 



(25) 



. _ (A,|Afe)2)] 

where — f3L/A and the numerical constants /„ (with 
n = 2,3,4) are given by the integrals 

(26) 







■ 4 
sm 77 



When we are interested in the structure of the supercon- 
ducting phase just below Tc(/), the comparison of the 
forth order terms in the free energy of different states 
should be made by calculating them at T = Tc{f). 
Thus, in the following we will need the values of In at 
77 = ric{l/2) — arccos(l/\/3), which are 

h = (15 -9a/27?c"^)/32« 0.0524, 

/3 = y3(7\/2r7c"^ - 9)/32 « 0.0737 , (27) 

h = (27- 13\/2?7c-^)/32 R:! 0.2424 . 
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B. A tripod of three wires 



C. The equal amplitude hypothesis 



For / = 1/2 the contribution to FnJ from a tripod 
formed by the three hnks (jak) [where ja = ia{k) with 
a = 1,2,3 are the three nearest neighbors of k numbered 
in the positive direction] after the substitution of Eq. IjlSf) 
can be rewritten as 



Fa 



Vl + V2 



3 

+r/2E (l^jJI^Jo+i|sinXj„j„+J^ - 

a=l 
3 



(28) 



where 



^3 



16 ^ 



V3 



whereas Xjj' the gauge-invariant phase differ- 

ences between the phases tpj of the order-parameter 
Aj = \Aj\exp(iipj) at neighboring six-fold coordinated 
sites. 



Xjj' = 'fir - Vj - 



-Xj' 



(29) 



It follows from Eq. p9|) that for all tripods, or, in other 
words, for all plaquettes of the triangular lattice formed 
by the six- fold coordinated sites. 



Let us now introduce the additional assumption (whose 
self-consistency is established in Appendix that the 
absolute values of the variables Aj are the same for all 
six-fold coordinated sites j, 

Aj — Aexp{iLpj) , (31) 

where A is real. In that case the contributions to 

J^it^=E^.^wH^) (32) 

k 

from the next but last term in Eq. (|28() coming from 
the neighboring tripods also cancel each other, and the 



expression for F^t) is reduced to 

FW=F4A4^K + n{x,.,„,J)] 

k 

where i^o = (9/2) i^i + 31^2, 



(33) 



3 \ ^ 3 

Viix}) = V{xi,X2,X3) = Vl I Vsinxa +;^2 Vsin^Xa 



a=l 



and 



E Xja]a+i = 7^/2 (mod 27r) 



(34) 



(35) 



a=l 



Since F^w is invariant with respect to the shift of any 
of the variables Xjj' by a multiple of 27r, here and below 
we for convenience assume that they all are reduced to 
the interval [— 7r,7r], in accordance with which the right 
hand side of Eq. (|35|l is written as 7r/2 (mod 27r) instead 
of — 37r/2, as it would follow from Eq. H30() . 

At the temperature of the superconducting transition 
in a fully frustrated wire network [that is, a,t rj = r]c = 

(1/2)] the values of the coefficients i^n are given by 



lyi = i(l + V2?7c"i)« 0.6201, 
V2 = 1(5 - 3^/2?7-l) « 0.0699, 

i/Q = 3 . 



(36) 



E 



Xjaja 



E^. 



Ja+l 



-37r/2 



(30) 



IV. MINIMIZATION OF THE FOURTH-ORDER 
CONTRIBUTION TO FREE ENERGY 



Since each plaquette of this lattice has a particular three- 
fold coordinated site in its center, the index k numbering 
such sites can be also used for numbering triangular pla- 
quettes. 

In the cyclic sums in Eq. (|28|) and analogous sums be- 
low j4 = ji. The last term in Eq. (|28|l can be omitted, 
since during summation over the whole lattice the two 
tripods adjacent to any link {jajb) always yield opposite 
contributions. 



A. A single triangle 

It is well known that the ground state of the 

antiferromagnetic XY model with triangular lattice can 
be found by minimizing the energy separately for each 
triangular plaquette and then matching these solutions 
with each other. 

For 1^1,1^2 > the minimum of V{xitX2iX3) oh an 
isolated triangle [under the constraint of the form l|35|l ] 



is achieved when two of the arguments of V{xi,X2,X3) 
are equal to each other, for example 



XI = X2 = -Ht) , X3 - 7T/2 + 2ij{t) , 



(37) 



another solution with the same value of V{xi, X2, Xa) be- 
ing 



XI = X2 = -[tt - , X3 = ^/2 - 2iP{t) 
where t = V2/V1 and 



(38) 



^(0 = arcsin^™±il^. (39) 
^ 4(1 + t) ^ ^ 

With increase of t from zero to infinity ipi^) continuously 
increases from arcsin [(VS — l)/2] ~ ti'/S to 7r/6. 

In an infinite system each variable Xjj' belongs si- 
multaneously to two triangles, but enters the function 
^(Xi:X2iX3) oil these two triangles with the opposite 
signs. Comparison of Eqs. (|37|) and H38|) with each other 

allows to conclude that it is impossible to minimize Fnv) 
by minimizing V{xi, X2, Xs) separately for each triangle. 



B. A periodic solution 

When variables Xjj' are reduced to a finite interval, the 
average value of J2l=i Xjaja+i should be equal to zero. 
This can be achieved if on one quarter of triangles the 
right hand side of Eq. (|35|l is equal to —3tt/2, whereas on 
all remaining triangles it is really equal to tt/2. Accord- 
ingly, the minimal supercell of a periodic set of variables 
Xjj' should consist of four triangles. 

The four-triangle supercell with the most symmetric 
(triangular) shape, but with the most general structure 
allowing for construction of a periodic state by a periodic 
repetition of this supercell, is shown in Fig. 2a. It can be 
described by the six variables Xi (with i = 1, . . . , 6) de- 
fined as shown in Fig. 2a and satisfying three independent 
constraints of the form 1351) . which can be chosen to be 



Xi +X2 +X3 = 7r/2 , 
-X3 + X4 + X5 = 7r/2 , 
-Xi - X5 + X6 = 7r/2 , 

the fourth constraint, 

- X2 - X4 - X6 = -37r/2 , 



(40) 
(41) 
(42) 



(43) 



following automatically from Eqs. (|4U|I - H42|I . 

The minimization of Fnw for this supercell with respect 
to the remaining three degrees of freedom shows that for 
Xi G [— TT, tt] the minimum is achieved when 



Xi = X4 = -7r/4 , X3 = , 
X2=X5=37r/4, X6=7r, 



(44) 



(a) 




(b) %i X. 



X 




Xi 



(c) 



AAA 

\AAA/ 
\/\A/ 



FIG. 2: Construction of periodic patterns minimizing the 
fourth-order contribution to free energy, (a) A possible choice 
for the most symmetric four triangle supercell. (b) An alter- 
native four triangle supercell. (c) The structure of the periodic 
solution obtained from free energy minimization with the su- 
percell shown in (a) or (b). Simple arrows correspond to phase 
differences Xjj' equal to ±7r/4, double arrows to ±37r/4, sim- 
ple lines to and wiggly lines to tt. 



or in one of the five other states which can be constructed 
from this state by permutations of the variables Xi- In 
these states on all triangles 



smxiaia 



(45) 



which means that on each triangle the first term of 
^({x}) reaches its absolute minimum (i.e., is equal to 
zero). Accordingly, the value of F,lw in these states does 
not depend on i^i. 



= 2(z/o + j^2)F4A47V . 



(46) 



Note that the supercell defined by Eqs. (|44|l consists 
of two pairs of equivalent (if one takes into account the 
equivalence of tt and — tt) triangles. Thus the actual size of 
the supercell has turned out to be two times smaller than 
it has been initially conjectured. But there was no way to 
predict this without really performing the minimization 
for the four-triangle supercell. 
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The same sohition (whose structure is shown in Fig 2c) 
can be also found by starting from the assumption that a 
periodic state is formed with the help of the four-triangle 
supercell with the different shape shown in Fig. 2b. In 
that case the constraints (I40II and (|41|1 retain their form, 
whereas in Eqs. H42|l and 143|l one should interchange xi 
and X4- For this supercell the minimum of F^w (for not 
too large ratio 1^2/1^1) is again achieved in the solution 
described by Eqs. (|44|) (or other equivalent solutions). 



It follows from Eq. that for \Aj 




a=l 



(47) 



Substitution of Eq. (|45ll into Eq. 1471) allows immediately 
to conclude that in the solution which we have found the 
absolute value of the order parameter on all three-fold 
coordinated sites has the same value as on the six-fold 
coordinated sites, 



|Afc| = |A,| = A. 



(48) 



In addition to gauge-invariant variables Xjj' defined on 
the bonds of triangular lattice, one can, naturally, also 
introduce the gauge-invariant phase differences defined 
on the bonds of the original dice lattice. 



jk ^ (fik - f] - Ajk , 



(49) 



where (pk = arg(Afc). As a consequence of Eq. 1)12(1 . the 
variables 9jk — ~0kj, which we assume to be reduced to 
the interval (— tt, tt), have to satisfy the constraint 



(50) 



When Eq. 1481) is fulfilled, the expression ||SJ) for the cur- 
rent in a link is reduced to 



where 



Ijk — Iq sin 9jk 
2e 7A2 



h 1^ sin 77 
Calculation of 

AjAlexp{iAjk) = A'^ exp{—i6jk 



(51) 



(52) 



(53) 



with the help of Eq. (|16|l demonstrates that in the con- 
sidered state the variables 9jk have the same three values 
(up to a permutation and a simultaneous change of sign) , 



9,k=±0i, ±82, Ti 



(54) 



on all tripods forming dice lattice. These values have to 
satisfy the constraints. 



leading to the automatic fulfillment of Eqs. (|50() . as well 
as the current conservation equation. 



L 09 = sin 9r, 



(56) 



which follows from Eq. H51() . As a consequence, they turn 
out to be exactly the same. 



n,3 



i(^1/%/3t1/%/6) , 6I2 = arccos(^l/%/3 



(57) 

as in the ground state of the fully frustrated XY model 
(FFXYM) with a dice lattice |2g, for which the current 



7r/4 , 6*1 + 6*3 = 7r/2 , 6*2 + 6*3 = 37r/4 , (55) one assumes |Aj 



conservation equation also has the form jSoJ. 

Thus, in terms ofOjk, the state which we have found (it 
is schematically shown in Fig. 3a) has exactly the same 
structure as one of the ground states of the FFXYM with 
a dice lattice. It has to be emphasized that the reasons for 
that are more subtle than a simple reduction of one model 
to the other. Firstly, in the case of a superconducting 
wire network the relation H48|) . the form of which seems 
to imply a possible reduction to XY model, is valid only 
in the minimum of free energy. Secondly, the substitution 
of Eq. igSJl into Eq. (|2SJ) gives 

4tre(^jfc) = ^l[2(/4 + I2) + 8/3 COSe^k + COs" Ojk] , 

(58) 

whereas in the case of the FFXYM in the expression for 
the energy of a bond the term proportional to cos^ 6 is 
simply absent, whereas the term proportional to cos 9 
has a coefficient of the opposite (negative) sign. Thus the 
two models do not become equivalent even if Eq. H48|) is 
artificially introduced as an additional assumption. 



C. Additional degeneracy 

The ground state of the FFXYM with a dice lattice is 
known to possess a well-developed accidental degeneracy, 
which can be described in terms of the formation of a net- 
work of zero-energy domain walls (26i | on the background 
of a periodic state. This construction can start from any 
of the four periodic states shown in Fig. 3 and allows to 
obtain, in particular, the three other states shown in that 
figure by introduction of such domain walls. 

For example, both state (c) and state (e) can be ob- 
tained from state (a) by introduction of a dense sequence 
of parallel zero-energy domain walls (of two different 
types). In that language state (g) can be described as the 
dense network of intersecting domain walls of two types. 
On the other hand, if one starts the construction from 
state (c), both state (a) and state (g) are formed by in- 
troduction of parallel domain walls, whereas the network 
of two types of walls corresponds to state (e) . 

The same set of states (described in a more detail in 
Ref. \2(^ minimizes for given A. In Appendix B we 
check that all these states are extremal not only when 
const, but also in the absence of this 
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(a) 




(e) 








(c) 




(g) 








FIG. 3: The four periodic states generating the class of degenerate states discussed in the text by adjunction of domain 
walls. The plaquettes with positive vorticities are shaded, and the three types of arrows correspond to the three possible 
gauge- invariant phase differences B-i, 62, S3 (modulo 2-k). Labels are chosen to match those introduced in Ref. 



constraint. The accidental degeneracy related to forma- 
tion of zero-energy domains walls gives rise to residual 
entropy proportional to the linear size of the system. 

In the framework of the description of different states 
in a network in terms of the gauge-invariant phase dif- 
ferences 9jk G (— tTjTt), all rhombic plaquettes a of a 
dice lattice can be considered as occupied by positive 
and negative half- vortices, whose vorticities nia, = ±1/2 
are given by 



±- 



1 



(59) 



In the family of states minimizing the half-vortices 
of the same sign always form triads with one "central" 
and two "edge" vortices [2^. The formation of domain 
walls which cost no free energy is related to the changes 
in the orientation and/or in the shape of these triads, but 
does not lead to formation of vortex clusters of any other 
size. 

At low temperatures {T <^ Tco) the free energy of a 
fully frustrated wire network with the dice lattice geom- 
etry (which then can be described in terms of the London 
approximation) is minimal for the same set of states, but 
with the slightly different values of 9a 

01 = 7r/12 , 62 = 7r/3 , 6*3 57r/12 , (60) 

satisfying, nonetheless, the same constraints (|55|l . 



D. Alternative solution 

The analysis of the supercell shown in Fig. 2b allows 

also to find a state which minimizes i^iw for large 1^2 ■ In 
the notation of Fig. 2b the structure of this state is given 
by 

Xi = X2 = X3 = -X4 = 7r/6; X5 = -Xe = 57r/6 . (61) 



It minimizes the value of the second term in Eq. I|34|) 

^(4) 1 • 

'ilw bemg 



separately for each triangle, the value of Fn 

3, 



FaA^N 



(62) 



This alternative state is characterized by even more de- 
veloped accidental degeneracy leading to extensive resid- 
ual entropy. Namely, the value of Fnw does not change 
if at an arbitrary number of sites j the variables ipj are 
shifted by tt. Note that this property holds not only at 
i^i = 0, when it trivially follows from T^({x}) being de- 
pendent only on cos 2^0, but also at finite values of vi. 

Comparison of Eq. (|62|l with Eq. (|46|l shows that 
the values of F^w in the two states become equal at 
1^2/1^1 = 3, whereas in our case, according to Eqs. 
1^2/1^1 « 0.1127 <C 3. At so low values of the ratio 1^2/1^1, 
the alternative solution discussed in this subsection is 
simply unstable. 
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V. MAGNETIC ENERGY 
A. An arbitrary network 

When currents in a two-dimensional wire network sat- 
isfy the current conservation equations, the current in 
each hnk can be expressed as a difference of so called 
mesh currents [i^, /™, associated with the plaquettes 
(meshes) of a network. Namely, the current in the link 
(jfc) can be written as the difference of the mesh currents 
associated with the two plaquettes (a and a') adjacent 
to this link. 



jk 



jra 



(63) 



The magnetic energy of the currents in the network, 
^'magn, Can be then expressed in terms of /™, 



(64) 



where a symmetric matrix L^p is usually called the mu- 
tual inductance matrix 40] . 

The diagonal elements of this matrix describe the self- 
inductances of current loops which can be associated 
with different plaquettes of the network and, accord- 
ingly, have to be positive. On the other hand, its non- 
diagonal elements describe the mutual inductances of 
non-intersecting coplanar current loops. Magnetic fields 
of such loops substract from each other and, therefore, 
the non-diagonal elements of L^fs have to be negative. In 
an infinite network the constraint 







(65) 



has to be satisfied. This ensures the invariance of -Emagn 
with respect to a possible redefinition of mesh currents, 
I™ /™ + (5/™ that leaves the physical currents in the 
links, Ijk, intact. For practical purposes it is convenient 
to define I™ in such a way that 







(66) 



The value of Lap depends only on the relative dispo- 
sition of the two plaquettes a and (3, and in the limit 
of infinitely thin wires can be found by calculating the 
double integral over their perimeters Fq, and F^, 



1 



(67) 



The expression for Lap given by Eq. H67|) in the case of 
a = (3 is logarithmically divergent, which means that 
Laa has always to be calculated more accurately, taking 
into account the finite width of the wires. The same is 
true for the value of Lap for neighboring plaquettes (hav- 
ing a common link). In the case of more distant neighbors 



(having only a common node or simply not touching each 
other) one can use Eq. (|67|l based on the assumption of in- 
finitely thin wires without encountering any divergences. 

Eq. (|67|l can be also rewritten as the double integral 
over the areas of the plaquettes a and /3. For a (3 



Lai: 



rp 



(68) 



which shows that \Lap\ rapidly decays with the growth of 
Rap, the distance between the centers of the plaquettes 
a and (3. For Rap ^ L 



Lap ~ — 



SaSp 
^aP 



(69) 



In any periodic state minimizing the free energy of a 
frustrated network, the plaquettes with negative and pos- 
itive values of /™ regularly alternate with each other, so 
-£'magn (normalized, for example, per a single plaquette) 
is given by a rapidly decaying sign alternating lattice 
sum. It allows one to expect that the main contribution 
to this sum comes from its largest terms, correspond- 
ing to the self-inductances of the plaquettes and the mu- 
tual inductances of rather close neighbors. Analogously, 
when comparing the magnetic energies of different degen- 
erate states minimizing the fourth order term in the free 
energy, the main contribution to the difference between 
them can be expected to come from the closest neighbors 
whose contributions do not cancel each other identically. 

Besides the proper energy of the magnetic field induced 
by currents and given by Eq. (|64|l . one also has to take 
into account the decrease of the superconducting free en- 
ergy related to the vector potential of this field. In the 
weak screening regime the sum of these two contribu- 
tions, Fmagn, differs from E^^gn only by sign, 



magn -^magn ; 

and, therefore, one has to maximize E'magn- 

B. Mutual inductances of dice lattice plaquettes 

Fig. 4 introduces the classification of neighbors for 
rhombic plaquettes of a dice lattice, which can be used for 
the natural reordering of summation in Eq. (|64|l . A cho- 
sen plaquette (which is shaded) has four nearest neigh- 
bors (denoted by 1), four next-to- nearest neighbors (de- 
noted by 2), two third neighbors (denoted by 3), etc., as 
shown in Fig. 4 up to sixth neighbors. In the following 
we denote the self-inductance of a plaquette Lq and the 
mutual inductance of a plaquette and its n-th neighbor 
(which is a negative quantity), L„ = — A„L. 

For n from 2 to 5 the calculation of L„ with the help 
of Eq. linijl or Eq. (jHUl gives 



A2 = 4V3-2%/7-2-151n3-41n(l-hV3) 
+ ln[8(4 + V7fi-1 + 2V7)(1 + 25\/7)2] 



(70) 
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A3 = 8\/7- 12V3 

+ ln[(27/4)(l + V3)iV(5 + 2V7f] , 
Xi ^ 12- 10\/3 + 2V7- (33/2)ln2- (97/4)ln3 (72) 

+ (1/2) ln[(l + V3)5°(4 - V7f\2 + ^7)"] , 
As 3 + 5V3-3V7- VT3- 121n3 (73) 

+ ln[(2 + V7){1 + 2%/7)(5 + 2V7f/{l + Vsf] 

+ (1/2) ln[(l + VT3){-2 + Vl3)\7 + ^13)72] , 



(71) acquires only two values, 



whereas numerically 

As ' 
A4 : 
As ^ 

Thus the decrease of |L„ 
n - 1. 



0.5569 
0.3637 
0.1671 
0.0723 



(74) 



with n is rather fast even for 



C. The form and the magnitude of magnetic energy 

In a general situation the values of mesh currents /™ 
should be found from Eqs. and, accordingly, are 

given by non-local linear combinations of link currents 
Ijk- Quite remarkably, substitution of Eq. (|16l) into Eq. 
^ allows to find that the values of mesh currents in a 
fully frustrated superconducting wire network with dice 
lattice geometry in the vicinity of Tc are given by the 
local expression. 



2e 7 
h ^ sin r] 



|A,||A,, 



V3 



(75) 



where I^, is the mesh current in the plaquette 



a,y = {jk'j'k"). 

In the family of states minimizing Fnw and described 
in Sec. IV the absolute value of 1"^, on all plaquettes 




FIG. 4: Classification of plaquettes according to their dis- 
tance from the shaded one. 



(76) 



where Iq is given by Eq. H52|) . The choice between I™ 
and 7™ is determined by whether the plaquette ajji is 
occupied by the central or edge vortex of a triad to which 
it belongs, whereas the sign of IJ^, by the sign of this 
vortex. The ratio g = I™ /I™ following from Eqs. I|76|l is 
equal to = 1/ V2. 

Substitution of Eq. (f75|l into Eq. shows that -Emagn 
is of the fourth order in Aj, and thus should be added 

(with the negative sign) to the forth-order term, Fnw, 
which has been minimized in Sec. IV. It follows from Eq. 
(|75|l that for |A^ | = A the main contribution to -Emagn 
[which is related to self-inductances of lattice plaquettes, 
Lq ^ 8Lln{L/d), where d is the thickness of the wires], 
can be written as 



magn 



(33') 



X33' 



(77) 



where E^^ — LoIq/Gc^- 

The V 
given by 



The value of the coefficient E^, which at T = Tc is 



E„ 



nc 



(78) 



should be compared with FjA"* = {L/4)l3A^. With the 
help of Eq. HA5|I and Eq. (|A6|) one obtains 



E„ 



F4A4 K^eiTc) ' 



(79) 



where k is the Ginzburg-Landau parameter of the ma- 
terial from which the wires are made, and a their cross- 
section area. Since the mean field phase transition in a 
fully frustrated network takes place when £,{Tc) ~ L, Eq. 
(|75)l can be rewritten as 



Fn 



F4A4 



(80) 



Thus, in order to get Fmagn ^ 7nw (the weak screen- 
ing regime) one should have 

a < K^L^ . 

In other terms the same condition can be rewritten as 

AeffK(r,) - L] » i , 

where Aeff(^) ^ {n^YL/a is the effective penetration 
depth for the magnetic field in a network. In such a case 
Fniagn Can bc treated as a small correction, which is rel- 
evant only for the removal of the accidental degeneracy 
between different states minimizing F^w . In the opposite 
limit finding the structure of the ordered state becomes 
an even more complicated problem because one has to 
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minimize F^w + F^w + E^^^^^ taking into account the 
dependence of all three terms on the magnetic fields in- 
duced by the currents. Therefore, below we always as- 
sume a K^L^, that is that the wires are thin enough. 

Note, however, that the applicability of the mean 
field approach requires the wires to be not too thin, 
a ^ k'^L^ / Aunw{T), see Appendix A. The two conditions 
on a are compatible provided L <C Auniv(T), which is 
readily satisfied in experiments, since for a temperature 
around one Kelvin Auniv(T') is of the order of one cen- 
timeter, see Eq. HA8|I . 

It is clear from the form of Eq. H77I) that the substrac- 

tion of -Emagn from Fn^ does not change the functional 
form of the fourth order term (expressed in the terms of 
the variables Xjj')j but leads only to a small increase of 
the coefficient 1^2 in Eq. H34|) . The next contribution to 
^magn; which is related to mutual inductances of neigh- 
boring plaquettes, has a more complicated structure. 



mag 11 



■2^2^cosxj„j„+i 

fe a=l 



+ l3a + 2 



(81) 



The non-equivalent contributions to magnetic energy 
of different periodic states can be characterized by the 
dimensionless parameter e defined by the relation 



i^magn = 42gn + ^i'ign + ^ ' (85) 

The values of e for basic periodic states can be then ex- 
pressed in terms of the coefficients A„ as 



2c/ X2 - (1 + 2.g2)A3 + (2 - 4g^)X^ 
+ 4g'A5 + . . . , 



(86) 
(87) 



= (l + ff')A2-ff^A3 + (l-35^)A4 

- (2 - 2.92)A5 + . . . , 
ee = -252A2 - (1 - 2.g2)A3 + 2A4 - VA5 + 
eg = (l-3.g')A2+.g2A3 + (l + .g2)A4 

- (2 - 252)A5 + . . . . (89) 

Substitution of the values of A„ given by Eqs. H74(l into 
Eqs. l|S^ - (|5^ shows that for \g\ < 1 the state (c) maxi- 
mizes the n-th approximation to e, 



However, in all the states minimizing Fnw , the three vari- 
ables Xjj' on each triangle are always given by 

Xi,2.3 = vr/2 T 7r/2, tt^/A, ±3^4 , (82) 

as a consequence of which the sum over a in Eq. H81|l is 
equal to —1/2 for all k and, therefore, the substraction 

of E'magn from Fnw — -Emagn also docs not lift the degen- 
eracy. To be sure of that we have checked also that all 
the states remain extremal with respect 

to variations of ipj even when one includes into analysis 

the extra terms obtained by the variation of i?magn. 



D. Selection of the state by magnetic energy 

In a more general situation, when the ratio g = I^/IJI^ 
is kept as a free parameter, the first two contributions to 
the magnetic energy can be written as 



magn 



magn 



Lo 
2 



' Tm 



{l + 2g^)N , 



2g'N 



(83) 
(84) 



and, naturally, are also the same for all the states mini- 
mizmg Piiw . 

(2) 

The difference appears when considering Fmagn, the 
contribution to Fmagn coming from the mutual induc- 
tances of the plaquettes which are next-to-nearest neigh- 

(2) 

bors of each other. Comparison of Fmagn for the four 
periodic states the structure of which is shown in Fig. 3 

(2) 

shows that for any g with I5I < 1 the minimum of E'^i 

gn 

is achieved in the state (e) and the maximum in the state 
(c). 



/ y ^magn 



(90) 



not only for n — 2 (as it follows from the previous para- 
graph), but also for all other values of n it has been pos- 
sible to check (up to n = 5). Since the value of A„ rapidly 
decreases with the increase of n, the same conclusion can 
be expected to be valid in the limit of n — > 00. Table 
1 illustrates the dependence of e*-"^ on n for four basic 
periodic states in the case of g = 1/ V2. 

In a Josephson junction array the value of the cur- 
rent in a junction is also given by Eq. I|51|l . where now 
Iq = {2e/h)J (J being the coupling constant of a single 
junction). As a consequence, the results of this section are 
applicable to a fully frustrated Josephson junction array 
as well. The value of the parameter g in the array is also 
equal to 51, because its ground state is characterized by 
exactly the same values of the variables 9a, Eq. l(S7|l. as 
in a fully frustrated wire network just below Tc. 

On the other hand, for T <C Tco, when Ijk oc 6*^^, 
the value of 17 in a fully frustrated wire network is equal 
to go = 3/5, as can be found from Eqs. It can be 

expected that with the decrease in temperature the value 
of (7 in a network continuously decreases from 171 to go- 
However, the conclusion on the selection of the state (c) 
by the magnetic energy is valid not only for go < g < gi, 
but in the whole interval — 1 < g < 1. 



VI. CONCLUSION 

The main result of this work is that the superconduct- 
ing state in a fully frustrated wire network with the dice 
lattice geometry exhibits the same set of degenerate spa- 
cial patterns minimizing the free energy in the two lim- 
iting cases when the temperature is either low compared 
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to the critical temperature Tc (London limit) or close to 
Tc (Ginzburg-Landau limit). This conclusion is quite in- 
teresting, given the fact that the system is described by 
two rather different models in these two limits. In the 
London limit, the amplitude fluctuations of the super- 
conducting order parameter are negligible, and the cor- 
responding fully frustrated XY model has been analyzed 
in 26]. In the vicinity of T^, we used the variational ap- 
proach pioneered by Abrikosov, where the spacial vari- 
ations of the complex order parameter are constrained 
to live in the subspace of unstable modes for the corre- 
sponding linearized Ginzburg-Landau equations. 

Our interest in this problem had been stimulated by 
the fact that for a fully frustrated network with the dice 
lattice geometry, this subspace has unusually high degen- 
eracy (as a consequence of the localized nature of modes) 
and includes a finite fraction (one third) of the total num- 
ber of modes Remarkably, non linear effects se- 
lect particular linear combinations of these spacially lo- 
calized states which reproduce precisely the current pat- 
terns obtained for the pure XY limit |2y|. 

In the second part of this article, we have compared 
magnetic energies of current patterns in different periodic 
states minimizing the Ginzburg-Landau functional. The 
dominant contribution to this degeneracy lifting mecha- 
nism is due to interactions of current loops which can be 
associated with second neighbor plaquettes. It favors the 
periodic pattern in which the triads of positive and neg- 
ative vortices have three different orientations [state (c) 
in Fig. The same conclusion is valid also in London 
limit and in the case of a Josephson junction array with 
the same geometry. It can be hoped that decoration ex- 
periments performed in more equilibrium conditions than 
those of Refs. |3 12^ may reveal such an ordering. 

This work leaves several open questions. First, it is 
important to know if all the states constructed here are 
stable with respect to local fluctuations in both order pa- 
rameter amplitudes and phases. We have checked numer- 
ically that this is indeed the case for the periodic states 
shown on Fig.|31 as well as for configurations with a single 
domain wall between two such ordered states. Unfortu- 
nately, a simple (analytical) stability proof holding for 
the complete class of degenerate states is not available 
now. 

Second, alternative degeneracy lifting mechanisms 
should be analyzed, as has been done for a wire network 



n 


a 


c 


e 


g 


2 


0.5569 


0.8354 


-0.5569 


-0.2785 


3 


-0.1704 


0.6536 


-0.5569 


-0.0966 


4 


-0.1704 


0.5700 


-0.2227 


0.1540 


5 


-0.0258 


0.4977 


-0.3673 


0.0817 



TABLE I: Comparison of dimensionless parameter e charac- 
terizing the magnetic energy of four basic periodic states in 
different orders of approximation. 



with kagome geometry by Park and Huse |20|| . These au- 
thors have found that the dominant perturbation to the 
idealized Ginzburg-Landau description of a wire network 
arises from the finite width of wires, which removes the 
degeneracy between the two opposite orientations of the 
supercurrents in a given loop. This mechanism can be 
interpreted in the terms of magnetic field redistribution 
between the network plaquettes (see discussion in Ref. 
I41L where it has been named "hidden incommensurabil- 
ity"), and, accordingly, is effective only when a network 
contains non-equivalent plaquettes. Therefore, in the case 
of a dice lattice (formed by identical rhombic plaquettes) 
it cannot play a prominent role. 

Another degeneracy lifting mechanism is related with 
the free energy of fluctuations around various free-energy 
minima and will be the subject of a separate report. It 
is likely to be the dominant one in the vicinity of Tc, 
but with decrease in temperature becomes less and less 
important in comparison with magnetic interactions of 
currents analyzed in this work. 
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APPENDIX A: CONDITION FOR THE 
APPLICABILITY OF THE MEAN FIELD 
APPROACH 

In two-dimensional superconductors the most impor- 
tant fluctuations are the fluctuations of the order param- 
eter phase = arg(A). When Tc is approached from be- 
low, the free energy of the phase fluctuations in a single 
wire, 



(Al) 



can be characterized by the effective coupling constant J, 
which in the absence of external magnetic field is given 
by 



J{T) 



7 



L/3 LI3e{T) 



(A2) 



The fluctuations are of no importance when 
J(r) ^ fc^T, which means that the critical region 
corresponds to 



em 



> 



7 



LpksTco 



(A3) 
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Since the maximal deviation of Tc{f) from Tco is achieved 
when £,[Tc{f)] ~ L the critical region can be con- 

sidered as sufficiently narrow if 



7 



(A4) 



When the mean free path of electron is much smaller 
than the thickness of a wire, the values of the coefficients 
entering the Ginzburg-Landau functional for the wire are 
determined simply by the values of the analogous coeffi- 
cients for the material from which the wire is fabricated. 
Substitution of 



(3 « cr/3bu 



7 ~ 0-7bulk : 



where a is the cross section of a wire, and 



_ 167r3 2 2 

Pbulk — ^2 7bulk''^ 
^0 



(A5) 



(A6) 



where k is the ratio of the penetration depth and the 
coherence length, into Eq. (|A4|I allows one to rewrite this 
condition as 



(Tco)ct 



(A7) 



where 



(2) 



site j. If this property is satisfied, it is then possible to en- 
force the equilibrium condition d[Fi 

with equal amplitudes \Aj \ sim 
J2j l^jP (with a negative coefficient below the critical 
temperature of the network). 



+ Fi^]/d\A,\^0 
is proportional to 



Derivation of i^iw 
gives 



X]fe^nw(fc) with respect to |Aj 



dFnw 



6 

4j/2 ^ sin^ 

b=l 



— 2vi sin 

6=1 



yijhjb+ 1 



(Bl) 

where (with 6 — 1, ... ,6) are the six neighbors of j 
numbered in positive direction (see Fig. [SJ, and it is as- 
sumed that = ji . The contribution from the third term 
in Eq. H28I) vanishes from Eq. (|Bip as a consequence of 
Eq. H45|l . which is valid for any tripod in any of the states 
discussed in section Hvl 

It turns out that the first sum in Eq. IjBip is equal to 
2 for any j in any of the considered degenerate states. 
On the other hand, summation of Eq. (|45|l over the six 
tripods containing the given site j allows one to conclude 
that the second sum in Eq. 1B1|) is always equal to zero. 
Thus, the expression in the right hand side of Eq. ljBl|) 
does not depend on j for any configuration with | A ^ | — 
const described in section ITVI 



Au„iv(r) 



2cm-K 

T 



(A8) 



is the expression for the universal value of the two- 
dimensional penetration length at the temperature of 
the Berezinskii-Kosterlitz-Thouless phase transition in a 
two-dimensional superconductor |42| . In the right-hand 
side of Eq. (jA8|) the temperature should be expressed in 
Kelvin. 

In aluminum wire networks fabricated with electron 
beam lithography H H i, El HI T^o ~ 1-2 K, and 

therefore Auniv(T'co) ~ 1.7cni,_whereas (according to the 
estimates of Park and Huse 0) k ~ 1. 




APPENDIX B: CONSISTENCY OF THE EQUAL 
AMPLITUDE HYPOTHESIS 



We shall now check that for the class of states described 
in section IIVI in which the amplitudes |Aj | are the 
same for all six- fold coordinated sites j, the derivatives 
dF^ /d\Aj\ are independent of the sixfold coordinated 



FIG. 5: The six neighbors of the site j contributing to the 
sums in Eq. UBIL 
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